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ABSTRACT 

The  tetrachoric  series  is  a  technique  for  evaluating  multivariate  normal 
probabilities  frequently  cited  in  the  statistical  literature.  In  this  paper 
we  have  examined  the  convergence  properties  of  the  tetrachoric  series  and  have 
established  the  following. 

Fpr  orthant  probabilities,  the  tetrachoric  series  converges  if 

/  fho  5" -  j 

]  p  -  - 1  <  l/(k-l),  1  i  <  j  <_  k,  where  p .  .  are  the  correlation  coefficients  of 

I  J  1  J 

a  k-variate  normal  distribution.  The  tetrachoric  series  for  orthant  proba¬ 
bilities  diverges  whenever  k  is  even  and  p..  >  1 / ( k- 1 )  or  k  is  odd  and 

j  J 

°i j  >  ,/(k-2)>  1  i1  <  j  ik  • 

Other  specific  results  concerning  the  convergence  or  divergence  of  this 
series  are  also  given. 

The  principal  point  is  that  the  assertion  that  the  tetrachoric  series 
converges  for  all  k  >  2  and  all  p.-  such  that  the  correlation  matrix  is 

J  J 

positive  definite  is  false. 
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SIGNIFICANCE  AND  EXPLANATION 


A  standard  technique  for  evaluating  multivariate  normal 
probabilities,  widely  cited  in  the  literature,  is  known  as  the 
tetrachoric  series. 

In  this  report,  it  is  shown  that  this  technique  is  defective 
in  that  the  series  will  sometimes  diverge. 


The  responsibility  for  the  wording  and  views  expressed  int'h'is  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


THE  USE  OF  THE  TETRACHORIC  SERIES 
FOR  EVALUATING  MULTIVARIATE  NORMAL  PROBABILITIES 


Bernard  Harris  and  Andrew  P,  Soms 
I.  INTRODUCTION  AND  SUMMARY 

Let  X  =  (Xj .Xg,. . . ,X^) ,  k  >  2,  be  a  normally  distributed  random  vector 
with  zero  means  and  unit  variances,  that  is,  X  has  the  probability  density 
function 

f~(xrx2 . xk)  =  (2ir)'k/2|RrV5sxTR  lx  ,  (1.1) 

X 

where  x  =  (Xj  ,x2». . . »xk)^,  A^  denotes  the  transpose  of  the  matrix  A,  and  R.the 
correlation, matrix  is  a  k  x  k  positive  definite  symmetric  matrix  with  elements 
p^.  and  =  1,  i  =  1,2 . k.  Further,  let 

♦(x)  =  (2tt)~^  e‘x?/2  , 

the  standard  normal  probability  density  function, and  let 

x 

<J>(x)  =  /  <fr(t)dt. 
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In  addition,  let 


=  ( 1  -<*>(  x)  )/4>(  x) 


and  for  i  >0,  let 


H.(x) 


=  ex  /2 


(-D1 


)• 


For  1  ^0,  H^x)  are  the  Hermite  polynomials. 
In  order  to  evaluate 


J|-  X  i  >h  i ,  X2>h  2  x  •  •  •  *  X  |^>h  }  =  / 

hk 


f~( X, ,Xp , 
X  1 


'■‘•G.  Kendall  [6]  proposed  writing 


P{X1>h] ,X2>h2,. . .  .Xk>hk;  = 


n,2=0 


l 


"y 

U 

n . . ! 


k 

n  h 


nk-1  k=0  l<i<j<k  ij‘  i=l  i 


(1.5) 


where  n.  =  Yn. .  +  Yn . . . 
i  L  ij  ^  ji 

j  •'  i < j  j:j<i 

The  right  hand  side  of  (1.5)  is  called  the  tetrachoric  series.  The 
particular  case  obtained  by  ordering  the  terms  of  (1.5)  with  respect  to  in¬ 
creasing  values  of  £n..,  namely 

"  J 


l  l 

ti=0  n 


I 


3.1J 

n .  . ! 


12 


Vl,k  U 


I  H  ,(hj»(hf), 
i=l  ni  1  1  1 


(1.6) 


where  n..  >0  and  £n..  =  N,  is  known  as  the  Kibble  series  or  the  tetrachoric 

*  J  *  J 

series  with  Kibble  ordering,  since  W.F.  Kibble  [9]  obtained  the  corres¬ 
ponding  series  for  the  multivariate  normal  density  function.  Note  that 
(1.6)  is  the  result  that  would  be  obtained  upon  integrating  the  series  given 
in  [9]  term-by-term.  Consequently,  to  avoid  ambiguity,  we  will  refer  to 
the  series  obtained  by  Kibble  in  [9]  as  the  Kibble  series  and  we  will  call 
the  corresponding  series  for  the  multivariate  normal  distribution  the 

-3- 


tetrachoric  series  with  Kibble  ordering. 


In  the  sequel,  to  simplify  notation,  the  right  hand  side  of  (1.5)  will 

be  denoted  by  Tk(h^  .h^,. . .h^)  and  Tk(0 ,0 . 0)  will  be  denoted  by  T^. 

Similarly,  the  right  hand  side  of  (1.6)  will  be  denoted  by  T£(h-j ^ > • • • .h^) 
and  T£.  Also,  sums  or  products  on  (i,j)  with  1  <  i  <  j  <  k  will  be  indexed 
by  i  <  j. 

The  tetrachoric  series  is  widely  employed  in  applications  of  multi¬ 
variate  analysis  and  is  quoted  in  many  standard  books  on  statistics  and  in 
many  papers  as  a  suggested  method  for  evaluating  (1.4),  often  with  the 
comment  that  convergence  may  be  slow  unless  |p..|  is  small  for  all  i,j, 

*  J 

i  f  j.  Specifically,  the  reader  is  referred  to  T.W.  Anderson  [1],  page  19, 
where  the  tetrachoric  series  is  suggested  as  a  way  of  evaluating  multi¬ 
variate  normal  probabilities.  In  R.E.  Barlow,  D.J.  Bartholomew,  J.M. 
Bremner,  and  H.D.  Brunk  [2],  page  137,  the  tetrachoric  series  is  mentioned 
as  one  of  several  possible  ways  to  evaluate  orthant  probabilities.  They 
comment  that  the  convergence  is  not  "fast  enough  for  practical  use  unless 
the  p..‘s  are  small".  R.L.  Plackett  [12]  comments  that  "although  the 

*  J 


tetrachoric  series  will  always  converge,  it  does  so  very  slowly  when  the 


absolute  values  of  the  correlation  coefficients  are  near  unity".  The  Kibble 

and  tetrachoric  series  are  described  in  connection  with  orthant  probabilities 

in  N.L.  Johnson  and  S.  Kotz  [5],  pages  44-46,  where  they  note  that  "these  series 

converge  very  slowly  unless  all  the  p..'s  are  small".  M.G.  Kendall  and  A. 

"i  J 

Stuart  [8]  describe  the  tetrachoric  series  on  pages  352-353  and  assert  that 

the  series  for  the  orthant  probabilities  always  converges,  but  only  slowly 

if  the  p..  are  not  small.  G.P.  Steck  [13]  described  the  Kendall  technique 
*  J 

noting  that  the  resulting  series  converges  slowly  when  the  p.^  are  large. 

S.S.  Gupta  [4]  discussed  both  the  results  of  Kendall  and  Kibble  for  the 
tri variate  normal  distribution,  remarking  that  the  series  converge  only  very 
slowly  for  high  values  of  |p..|.  In  [7],  M.G.  Kendall  employed  the  tetra- 

*  J 

choric  series  for  orthant  probabilities  to  study  the  distribution  of  upruns 
(sequences  of  increasing  observations)  in  a  time  series.  For  this  purpose, 
he  specifically  gives  the  terms  of  the  tetrachoric  series  for  the  orthant 
probabilities  for  k  =  2,3,4.  He  says  that  "the  expressions  are  not  as 
difficult  as  they  look"  and  that  "they  converge  fairly  quickly  for  damped 
autoregressive  series".  Kendall  says  further  that  "they  are,  he  thinks, 


amenable  to  calculation". 


In  this  paper,  we  have  reexamined  the  convergence  of  the  tetrachoric 


series  and  have  ascertained  that  under  specified  conditions  on  the  corre¬ 
lation  matrix  R,  the  tetrachoric  series  will  in  fact  diverge. 

The  development  in  this  paper  is  as  follows.  In  section  two,  we  dis¬ 
cuss  the  Kibble  series  for  the  multivariate  normal  density  function,  which 
is  an  intermediate  step  in  arriving  at  the  tetrachoric  series.  This  will 
permit  us  to  subsequently  indicate  the  difficulty  with  the  tetrachoric 
series.  A  sufficient  condition  for  the  convergence  of  the  tetrachoric  series 
is  given  in  section  three.  The  third  section  is  also  devoted  to  illustrations 
of  the  non-convergence  of  the  tetrachoric  series. 

2-  The  Kibble  Series  for  the  Multivariate 
Normal  Density  Function 

F.G.  Mehler  [10]  derived  series  expansion  for  the  density  function 
of  the  bivariate  normal  distribution.  This  was  subsequently  extended  to  the 
general  k-variate  density  (1.1)  by  W.F.  Kibble  [9].  Unfortunately,  Kibble's 
Kaper  contains  some  defects  and  as  we  will  show,  the  Kibble  series  does  not 
necessarily  converge.  Since  this  is  an  intermediate  step  in  arriving  at 


the  tetrachoric  series,  it  is  desirable  to  rederive  the  Kibble  series,  which 


will  permit  us  to  obtain  both  a  correct  result  as  well  as  to  indicate  the 
error  in  Kibble's  paper.  Accordingly,  we  first  establish  the  following  theorem. 

Theorem  i •  For  -i/(k-i)  < < i/(k-i) ,  l  <  i  <  j  <  k. 


f~(x,  ,Xp , . . . ,x. )  £ 

X  n!2=0 


ij 


k-1  ,k 


=  0 


n 

i<j 


JLL 

’id 


n  h  (x.)<j>(x.) . 
i=l  ni  1  1 


(2.1) 


This  series  is  known  as  the  Kibble  series  and  converges  absolutely  if  the 
p.j j  satisfy  the  above  restrictions. 

Proof.  Let  ^(t]  ,t2, . . .  ,tR)  =  e~ht  Rt,  t  =  (t] ,t2, . . . ,tk)T. 
t-|  ,tg » -  - . , t^ )  is  the  characteristic  function  of  the  k-variate  normal  dis¬ 
tribution.  Then,  from  the  inversion  theorem. 


~(x 

X 


(27  y)'k  / 


eU  x  ^(trt2 


. .  ,tk)dt-j dt2- - -dtk . 


(2.2) 


Writing 


(t)  = 


-ht^  t- 


l  P 
i<i 


ij 


t.t.  -ht't 
i  j  -  e 


I 

N=0 


l--1 )_  (  y  p  t  t  )N 


(2.3) 


and  expanding  the  last  factor  in  a  multinomial  series,  we  obtain 


\p(t)  =  e"  2 


-ht  t 


n12  Vl,knl  "k 

Pi  9  •  •  u  ti  •  •  • 


y  (-nN  y  H2  ‘Mk-1  ,k  L1 
=n  n12 ! * - -nk-l  ,k  ! 


N=0 


(2.4) 


where  the  inner  sum  is  over  n..  >0  with  In.  .  =  N.  A  rearrangement  of 

*  J  •  J 

terms  yields 


*(t)  =  e'^ 


CO  CO 


I  l 


"12'° 


n 


(-p^i 


nij 


k  n; 

n  t, 


Ji  n  I  S' 

Vi,k=0i<J  1J'  i=1 


(2.5) 


1 


obviously  dominates  the  partial  sums  in  (2.4)  and  (2.5)  and  as  the  follow¬ 


ing  argument  shows,  (2.7)  is  integrable. 

Note  that 

OO  Co  ^ 

(2tt)  /  ...  /  e~2t  tg(t)dt1dt2  . .  .dt.  =  2k  |  E|  2P(  >0  , . . . , 

—  00  -OO 


where  (X-j  ,X2 , . . .  ,X^)  has  the  multivariate  normal  distribution  with 
Z  ^  =  (a1''),  =  1,  oij  =  -p,  i  ^  j  and  the  hypothesis  |p|  <  l/(k 
sures  that  Z  is  positive  definite.  Thus,  we  can  substitute  (2.4) 
into  (2.2)  and  interchange  the  order  of  integration  and  summation. 


f~(x1,x2,...,xk)  =  (2iT)'k  l 
X  n12=® 


I 


n 

i<j  nij! 


k 

n 

i=l 


“  -1t.x.-*st.  n. 

fe  1  1  Vdti 

— oo 


Using  (1.3),  the  integral  in  (2.9)  is  easily  evaluated  obtaining 


xk>°}, 

(2.8) 

-1)  in- 
or  (2.5) 
Thus 


(2.9) 
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•  Hn  (xi)4>(xi). 
i 


k 

The  first  part  of  the  conclusion  follows  upon  observing  that  £  n.  =  2  £  n--. 

V  i=l  i<J 

k  n.  ,\"ij 

Hence  n  (-1)  1  =  (-1)  <J  ,  verifying  (2.1).  The  absolute  convergence  follows 
i  =  l 

from  (2.7),  (2.8)  and  the  observation  that  convergence  for  | pi ^ |  <  l/(k-l) 

implies  absolute  convergence,  since  (2.1)  is  a  power  series. 

REMARKS •  In  [9],  Kibble  actually  ordered  the  terms  as  indicated  in  (1.6). 

This  is  equivalent  to  (2.1)  whenever  the  series  converges  absolutely. 

Kibble  did  not  assume  |p..|  <  1 / ( k- 1 )  and  merely  asserted  that  term-by-term 

^  d 

integration  is  permissible  and  that  the  series  would  then  be  absolutely  con¬ 
vergent  for  all  values  of  the  variables  if  R  is  positive  definite. 

We  now  reexamine  the  hypothesis  Ip^jI  <  1 / ( k- 1 )  to  show  that  if  this 
condition  is  violated,  the  Kibble  series  may  in  fact  diverge. 

If  k  =  2,  f~(x, ,x?)  is  the  bivariate  normal  probability  density 
X 


-10- 


function  and  the  condition  of  Theorem  1  always  holds.  Theorem  1  is  pre¬ 
cisely  Mehler's  theorem  in  that  case.  Hence  consider  k  >  2. 

For  n  >  0, 


H  (x)  =  n! 
nv  ' 


r-i 

|  (-i)V-2" 

m=0  2mm! (n-2m) ! 


(2.10) 


so  that 


H„(0) 


r 

\ 


0 

n!(-l)n/2 
2n/2(n/2) ! 


V 


n  odd 

n  even  . 


(2.11) 


We  will  now  calculate  f ~ ( 0 ,0, . . . ,0) ,  using  the  Kibble  series,  for  the 

X 

special  case  p.  .  =  p,  i  f  j,  and  compare  this  with  (1.1),  which  gives 

*  *J 

f~(0 ,0 , . . .  ,0)  =  (27T)-k/2|R|--  . 

X 

Substituting  (2.11)  into  (2.1)  we  get 


M0,0,...,0)  =  (27r)‘k/2  °i 
X  n12=^ 


l 

nk-l  ,k 


k 

n  - 


n . ! 

_  _ 

iTwimf  n. .! 
i<j  J 


(-P/2) 


k 

l  n./2 

i=l  1 
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(2.12) 


Thus  for  any  pQ  for  which  (2.12)  converges,  it  will  converge  absolutely  for 
any  p  with  |p|  <  pQ.  For  such  p,  we  can  rearrange  the  order  of  terms  in 
(2.12)  so  that  it  is  a  power  series.  Further,  if  p..  =  p,  i  f  j,  it  is 

*  J 

easily  seen  that 


I R |  =  (l+(k-l)p)(l-p)k'1. 


-k/2 

Hence,  the  power  series  expansion  for  (2tt)  |R|  must  be  the  same  as 

the  power  series  indicated  by  (2.12).  However,  the  power  series  for  |  R  ^  has  a 
radius  of  convergence  of  l/(k-l)  which  implies  the  divergence  of  (2.1)  when¬ 
ever  | p. . |  >  l/(k-l),  1  <  i  <  j  <  k.  Thus,  we  have  shown  that  in  the 

*  J 

equicorrelated  case,  Theorem  1  can  not  be  improved  upon. 


3.  The  Convergence  of  the  Tetrachoric  Series 

Kendall  [6]  obtained  the  tetrachoric  series  by  integrating  the 
characteristic  function  term-by-term  as  in  (2.9).  Then  he  integrated  the  re¬ 
sulting  expression  (2.1)  term-by-term  obtaining 


-12- 


CD 


CD 


(3.1) 


/  •••  /  f ~( x-j , X£ , . , .  ,x^ )  dx,  dx„ . .  .dx^  Tj,(h,  .h^ , . . .  ,h^) . 

h,  h,  X 

k  I 

We  now  discuss  the  convergence  of  the  tetrachoric  series. 

THEOREM  2.  Let  h^.h^.-.-h^  be  arbitrary  real  numbers.  Then  if 
|p.j |  <  l/(k-l ) ,  1  <  i  <  j  <  k,  the  tetrachoric  series  Tk(h^  .h^, . . . ,hk) 
converges  absolutely. 

Proof.  From  Erd£lyi,  Magnus,  Oberhettinger,  Tricomi  [3],  p.  208, 
for  n  >  0 


|Hn(x)|  <  ce 


(3.2) 


where  c  ~  1.086.  Then,  from  (2.1), 


Let  p  =  max  J p .  . j  <  l/(k-l).  Then, 
i<j  J 


Ip,/"  k 

"--nV- 

1<J  1J  1  = 


ip^rj  « 


n  (n.!)2  =  (  n  r-  n  j~  n  n.s) 
i=1  i<j  ni j ‘  i>j  ji'  i=l 


n  n.!  n  — -i — n-  n  — ^ — r 
i=1  j:i<J  ij'  j : i > j  nji' 


The  expression  inside  the  bracket  is  a  term  of  the  multinomial  expansion  of 

(  l  I  Pi  4 1  +  I  I  P-s  ,*  I )  1  -  Thus 
j:i<j  J  j : i > j  J 


I P,-  I  1J  k  ^  En./2  / 

l  |p,,l  +  l  Ip„i>  ’  idk-np)1' 

i<J  ij  i=l  j :  i<j  J  j :  i  >  j 


since  l  n  -  2  y  n...  Substituting  this  into  (3.3),  we  have 


c  e 


(2ir)k/2  i<j  ni ^=0 


n  l  ((k-i)p) 


ij 


-  I  x  /4 
c.ke  1  =  1 


(2ir) 


YJ2~  O-(k-l)p) 


-\-k(k-l)/2 


(3.4) 


Since  (3.4)  dominates  all  partial  sums  of  the  Kibble  series  and  is  ob¬ 
viously  integrable,  the  dominated  convergence  theorem  applies  and  the  Kibble 
series  is  integrable  for  |p..|  <  l/(k-l).  Hence, 

*  J 


P{X1>h1,...,x|f>hk}  = 


co 


00 


l  ...  I  IT 
"12=0  "k-1 ,k=0 


k  °° 

n  /  H  (x-)4>(x.) dx . . 
i=l  h.  ni  1  11 


For  n^  >  0,  from  (1.3) 


/  Hn_(x)$(x)dx  =  Hn._i(hi)4>(hi) 


(3.5) 


and  for  =  0,  since  Hq(x)  =  1 ,  we  get 


7  Ho(x)<}>(x)dx  =  /  tj>(x)dx  =  l-$(h)  =  4>(h) » 


as  required  by  (1.2),  and  (3.1)  follows. 

To  show  the  absolute  convergence  of  Tk(h] >h2>  •  • -hk) »  we  first  note  that 
for  n  =  0 


|H_i (*) |  <  3e 


x2/2 


and  for  n  >  0 


|Hn_-,(x)|  <  cex  /4((n-l)!)H  <  3( nj ) 


h/ll 


Consequently,  by  a  minor  modification  of  the  argument  leading  to  (3.4) 


|Tk(hrh2,...,hk)|  1— |k/2  (1-(k-l)p)~k(k~1)/2> 


-16- 


for  |p.jj|  <l/(k-l) ,  establishing  the  absolute  convergence. 

We  now  investigate  the  possible  relaxation  of  the  hypothesis 

|  P  -  - 1  <  1/(k-l),  1  £  i  <  j  <  k.  For  this  purpose,  we  first  investigate  the 
*  \) 

computation  of  the  orthant  probabilities  for  k  =  3  using  the  Kibble  ordering. 
Thus,  we  compute 


(2’r3/U  1  V"23‘’ 


n^  ni3  n23 

P12  P1 3  p23 


(0), 


From  (2.11),  in  order  that  the  product  H  H  (0)  not  vanish,  we  must  have 

i=l  ni-l 

every  n.  =  0  or  an  odd  positive  integer.  We  show  that,  for  k  =  3,  the  only 
non-vanishing  terms  are  n^  =  n-|3  =  n23  =  0  or  two  of  ni 2 ,nl 3 ’n23  are  zero 
and  the  remaining  one  is  odd.  To  see  this,  note  that  if  one  of  n-j 2 *nl 3 *n23 
is  zero,  say  n^2>  then  n^  and  n23  are  odd  and  n^  +  n23  is  even  and  the 
term  vanishes.  If  none  are  zero  and  n,2  +  n^3  is  odd,  then  either  n^2 
or  n^3  must  be  even.  Assume  n^2  is  even,  then  n23  is  odd  and  n13  +  n23 
is  even,  resulting  in  a  vanishing  term.  Consequently,  for  k  =  3,  the  non¬ 
vanishing  terms  are  (0,0,0),  (1,0,0),  (3,0,0),...  and  their  permutations. 
Therefore 
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r 


T0  =  (h)*  +  (  2tt  ) 


-3/2 


l 

m=0 


2m+l ,  2m+l^  2rrH-l 


J]2 


+P1 3  +p 


23 


(2m+l ) ! 


H2m(°)H-l 


(0) 


4tt 


I 

m=0 


2m+l _L_2m+l ,  2m+l 


12 


+P1 3  +p 


23 


2m+l 


(2m)! 

22m(m! )2  _ 


■|  (cos'1p-|2+cos'lpi3+cos  lp23^ 

2  4tt 


the  well-known  formula  for  the  trivariate  orthant  probability;  also  the 
power  series  for  cos’^x  converges  for  |x|  <1.  Hence,  despite  the  di¬ 
vergence  of  the  series  for  the  density  function  at  0  for  p  >  h,  the 
tetrachoric  series  for  the  trivariate  orthant  probabilities  converges  for 
all  pf .. 

The  analysis  of  the  convergence  of  the  tetrachoric  series  for  orthant 
probabilities  for  k  >  3  is  somewhat  more  complicated.  To  aid  in  this  analysi 
we  introduce  the  following  lemma. 

I_EMMA  If  k  >^2  is  even,  then  the  non-zero  terms  in  (1.5),  when 
hi  =0,  i  =  1 ,2, . . . ,k, are  those  terms  in  which  an  even  number  (including 
zero)  of  the  n^’s  are  zero. 
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If  k  >_  3  is  odd,  a  term  in  (1.5)  is  non-zero,  when  =  0, 
i  =  l,2,...,k,  if  ana  only  if  an  odd  number  of  the  n^'s  are  zero. 

Proof.  From  (2.11),  (1.5)  and  (1.2),  a  term  in  the  tetrachoric  series 
for  the  orthant  probabilities  is  non-zero  if  and  only  if  each  n^  is  either 
zero  or  a  positive  odd  integer. 

Considering  such  terms,  if  k  >_  2  is  even  and  an  odd  number  of  the  n^'s 
k 

are  zero,  then  £  n.  =  2zn.  .  is  even  and  is  the  sum  of  an  odd  number  of 
i=l  1 

positive  odd  integers,  which  is  a  contradiction. 

Similarly,  if  k  >_  3  is  odd  and  an  even  number  of  the  n^ 's  are  zero, 
k 

then  V  n.  is  the  sum  of  an  odd  number  of  odd  numbers  and  can  not  be  even. 
i  =  l  1 

We  now  investigate  the  convergence  of  the  tetrachoric  series  for  k  >  3. 
LeMA  2.  Let  k  >  4  be  an  even  integer  and  let 

n12  =  n34  =  •••  =  nk-l  k  =  »  n-jj  =  2m  otherwise,  m  =  0,1,2,...,  . 

Then  if  p  >  1 /( k- 1 ) , 

N  ,  k  \  . 

k(p)  =  ~ 6 -  11  Hn  I  (0)J  (4>(0)k  (3.6) 

m>K  .n.n.  . !  M»1  ni  1  ' 

1<J  1J 
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tends  to  infinity  as  m  ■*  where  N  =  £  n.  . . 

i<j  J 

Proof.  Since  N  =  (2m(k-l)+l), 


§{2m(k-l)+l) 

M">k(P)  '  ((Ll)!)k/2((2»0!)k(k-2|7?  <H2m(k-1) 


(0))k(4»(0))k. 


Write 


|{2m(k-l)+l) 


Mm.k(p)  =  (2mt|)k/2(2m),k(k-1)/2  <H2m(k-l) 


(0))k(*(0))k. 


Then,  using  (2.11)  and  the  elementary  inequality 


(2*)he~nnri+h  <  n!  <  2(2tt)  W**5  ’ 


(3.7) 


we  get 
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which  clearly  tends  to  infinity,  whenever  p  >  l/(k-l).  This  leads  to  the 
following  therorem. 

Theorem  3.  The  tetrachoric  series  (1.5)  with  iu  =  0,  i  =  l,2,...,k, 
diverges  for  k  an  even  integer  whenever  >  p  >  l/(k-l),  1  <_  i  <  j  £  k. 

Proof.  The  proof  is  immediate  upon  observing  that  for  each  fixed  m, 

Mm  k^  a  term  in  the  tetrachoric  series  for  orthant  probabilities  when 

P-jj  =  p,  1  £  i  <  j  <  k.  Thus  we  have  exhibited  a  divergent  subsequence  of 
terms,  which  clearly  precludes  convergence.  The  same  lower  bound  (3.8) 
applies  if  p^  >  p,  1  £  i  <  j  £  k,  yielding  the  same  conclusion  in  this  case. 

We  now  consider  the  tetrachoric  series  with  Kibble  ordering  for  the 
orthant  probabilities  and  even  values  of  k,  restricting  to  the  case 
P-jj  =  P.  1  1  i  <  j  £  k.  This  is  a  power  series  in  p,  of  the  form 
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v  N 

l  i.P  and  the  previous  argument  does  not  show  that  a,.  ,  does  not  tend 
N=0  ",K 


to  zero,  since 


“n  k  =  l  n  (nil-0_,(«(°))  n  H  ,(0)  (3.9) 

nj j>0  i<j  1J  i=l  V 

Vn.  .=N 
L  1J 

and  in  Lemma  2  we  made  a  specific  choice  of  the  n^.'s,  rather  than  computing 

J 

“N.k- 

Theorem  4.  The  tetrachoric  series  with  Kibble  ordering  (1.6)  for  orthant 
probabilities  diverges  for  even  k  >  4,  whenever  p,. .  ^  p  >  l/(k-l), 

1  <_  i  <  j  5  k. 

Proof.  Let  N  =  ^-(2m(k-l  )  +  l ) .  The  Mm  ^ ( p)  (3.6)  is  a  term  in  (3.9) 

‘or  that  N.  From  Lemma  1,  the  non-zero  terms  of  (3.9)  all  have  an  even 
number  of  n.'s  equal  to  zero. 

If  all  n.j  are  odd,  then  from  (2.11),  the  sign  of  such  terms  is  given  by 


(ln.-k)/2 
(-1)  1 


(-1) 


N-k/2 
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and  depends  only  on  N.  Thus,  for  fixed  N,  all  terms  with  the  opposite 


sign  must  be  among  those  terms  with  a  positive  even  number  of  n.'s  equal 
to  zero.  Let  ^  be  the  sum  of  all  terms  in  (3.9)  with  no  n^'s  zero  and 
let  ^  be  the  sum  of  all  terms  in  (3.9)  with  at  least  two  n^'s  equal  to 
zero.  Then,  for  p  >  0, 


0  iMm,k(p)  i^,k(p)pN’  (3.10) 

where  pp)  >  0,  since  N  =  |f2m(k-l)+l]  and  M  -  iy  is  necessarily  even. 
Obviously,  o„ik  =  yHik  ♦  SNJ(. 

Consider  the  tetrachoric  series  with  Kibble  ordering  for  £  =  k  -  2. 

Since  each  term  in  ^  has  at  least  two  =  0,  each  such  term  appears  in 
the  tetrachoric  series  with  Kibble  ordering  for  £  variables  with  the 

p 

omission  of  the  factors  (H  ^(0)<J>(0))  .  Thus, 


oo  oo 

I  |6N  kl (p|N  <  ❖2(0)H^1(0)(^)  l  . 
N=0  5  n-^O 


Vl,£=0 


*£(o)  n 


JpI 


ni  j 


n  T~  n  |H  ^(0) 
i<j  °i j '  1*1  ni  1 


N 

Since  for  l/(k-l)  <  p  <  l/(k-3),  kP  converges  absolutely  by  Theorem  2, 

N 

3^  kP  +  0  as  N  +  ®  and  Mm  ^(p)  -*■  °°.  Further,  (3.10)  implies  that 

°°  N 

l  cu,  .p™  diverges  for  l/(k-l)<  p  <  l/(k-3),  but  this  is  a  power  series 
N=0  * 

and  hence  divergence  is  established  for  all  p  >  l/(k-l).  The  conclusion 
now  follows  trivially. 

There  is  an  extensive  literature  on  the  evaluation  of  orthant  pro¬ 
babilities  for  k  =  4.  Many  references  may  be  found  in  N.L.  Johnson  and 
S.  Kotz  [5],  pages  53-58.  One  of  these  methods  merits  coiment  here  since  it 
specifically  employs  the  tetrachoric  series  with  Kibble  ordering.  In 
P.A.P.  Moran  [11],  the  tetravariate  orthant  probability  is  used  to  calculate 
the  variance  of  the  Spearman  rank  correlation  coefficient.  To  obtain  the 
orthant  probability,  Moran  rederived  the  tetrachoric  series  and  explicitly 
calculated  the  first  few  terms  of  the  tetrachoric  series  with  Kibble  ordering. 
He  made  no  comment  about  the  validity  of  the  series,  but  did  comment  that 
for  the  problem  at  hand  calculation  is  tedious. 

The  convergence  properties  of  the  tetrachoric  series  for  orthant  pro¬ 
babilities  appear  to  be  more  difficult  to  investigate  when  k  is  odd.  These 
properties  are  discussed  in  the  following  theorems. 
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THEOREM  Let  k  be  an  odd  integer^  3  and  order  all  the  correla¬ 
tion  coefficients  by  absolute  magnitude,  that  is, |p^|  >_  } |  >  ...  >  |  p  :, 
q  =  so  that  has  the  vth  largest  magnitude  among  the  p.^ ,  1  <  i  <  j  <  k. 
If  the  tetrachoric  series  (tetrachoric  series  with  Kibble  ordering)  for 
orthant  probabilities  converges  absolutely  for  l  =  k-1  with  correlation  co¬ 
efficients  pv,  v  =  1,2, ...,(2),  then  the  tetrachoric  series  (tetrachoric 
series  with  Kibble  ordering)  for  the  orthant  probabilities  converges 
absolutely  for  k. 

Proof.  We  have 


T 


k 


<  l  •••  l  M0))k  n 

n12’°  nk-l,k=0  i<j 


n . . 
I  U 


ij 


n . 


U 


k 

n 

i=l 


(3.11} 


From  Lemma  1,  every  non-zero  term  in  must  have  at  least  one  n^  equal  to 


zero.  Thus,  successively  setting  each  n^  equal  to  zero,  we  get 


M'j 


T.  1 H  i(°)  l  l  (4>(o))k  n 
v=l  n .  .>0  i<j 


ij 


-  "  I”  .,(0)1  i 
"«■  1=1  1 


1 


(3.12) 


the  right  side  of  (3.12)  dominates  (3.11)  since  we  have  enumerated  terms 
with  more  than  one  n..  =  0  more  than  once.  For  each  fixed  v,  the  inner  sum  dom¬ 
inates  <j>(0)Tk  -] .  Thus  (3.11)  converges  absolutely  whenever  each  of  the 
k  series  in  (3.12)  converges  absolutely.  Each  such  series  depends  on 

(2)  correlation  coefficients.  For  each  series  in  (3.12)  if  is  not  in 

0 

the  set  of  (^)  correlation  coefficients  of  largest  magnitude,  we  replace  it 

by  one  that  does  not  appear  in  the  series  in  such  a  way  that  each  p  , 

Jl  - 

v  =  1,2, ...,(2)  appears  once.  Thus, 


ij 


Tk  <  k(4>(0))KH_1(0) 


n]2=0 


n 

:0  i<j 


n. 


U 


SL 

n 

i=l 


which  converges  by  hypothesis.  The  conclusion  for  the  tetrachoric  series 
s.ith  Kibble  ordering  follows  upon  rearrangement  of  the  terms. 

Remark.  We  have  actually  proved  the  absolute  convergence  of  the 
tetrachoric  series  whenever  each  of  the  k  series  in  (3.12)  converges 

absol utely. 

COROLLARY  l.  If  each  =  p,  then  the  tetrachoric  series  (tetrachoric 
series  with  Kibble  ordering)  for  the  orthant  probabilities  converges 
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absolutely  whenever  -l/(k-l)  <  p  <  l/(k-2),  for  k  an  odd  integer  5. 

Proof.  The  conclusion  follows  from  combining  Theorem  5  and  Theorem  2. 
COROLLARY  2.  If  k  =  3,  the  tetrachoric  series  (tetrachoric  series  with 
Kibble  ordering)  converges  absolutely. 

Proof.  This  is  immediate  from  Theorems  5  and  2. 

The  observations  utilized  in  Theorem  5  also  enable  us  to  establish  the 
following. 

THEOREM  6,  If  k  >  5  is  odd  and  if  p. .  >  1 / ( k -2 )  for  1  <  i  <  j  <  k, 

* 

then  the  tetrachoric  series  (tetrachoric  series  with  Kibble  ordering)  for 
the  orthant  probabilities  diverges. 

Proof.  Consider  the  tetrachoric  series  with  Kibble  ordering.  From 
Lemma  1,  the  contribution  to  k  (3.9)  is  zero  unless  an  odd  number  of  the 
n.'s  are  zero  and  the  remainder  are  positive  odd  integers.  Write 
ClN  k  =  Yn  k  +  k’  where  yn  k  are  these  terms  with  exactly  one  n^  =  0.  Then 

yn  k  =  k  l  n  (n  ..!)  1  ^(0)  H  (0)  n  H  ,(0), 
n .  .>0  i<j  1J  i=l  "i  1 

■  J 

Jnij=N 

M 

where  i  =  k  -  1  and  1  <  i  <  j  <  l.  From  Theorem  4,  IyN  diverges  for 
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p  >  1  / ( k-2 )  =  l/(£-l).  3^  corresponds  to  those  terms  with  at  least  3  or 

more  ni  =  0.  Each  such  term  occurs  in  the  tetrachoric  series  with  Kibble 
ordering  for  k  -  3  variables.  For  l/(k-2)  <  p  <  1 / ( k-4) ,  that  series  is 
absolutely  convergent.  Thus,  the  methodology  of  Theorem  4  applies  and  we 
obtain  the  conclusion. 

So  far  we  have  only  considered  the  case  of  orthant  probabilities. 
Therefore,  we  now  show  that  similar  results  obtain  in  a  somewhat  more  general 
case.  For  this  purpose,  we  require  the  following  lemma. 

Lemma  3.  Let  x  sir/2(k-l),  where  s  and  k  are  integers,  and 
k  >  2.  Then 


cos  ujx)  =  cos  15  x  .  H.)  ; 


does  not  tend  to  zero  as  m  tends  to  infinity. 

Proof.  Let  m.  =  4j2(k-l),  j  =  1,2 .  Then 

J 
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cos  u 


(x)  =  cos  [4j2(k-l)2-  x  -  2j2(k-l)2ir+p 

■l  % 

=  -  sin  [2j(k-l)(l  — ~ - o)  x]  . 

8j  (k-lT 


Since 


2j(k-'l)  "  43T^TI  <  2j(k‘1)(1‘  8j2-(j.1)Z)  2  <  2j(k‘1)’ 


for  any  e  >  0,  there  is  a  jQ  such  that  for  all  j  >  jQ 


|sin[2j(k-l)(l - 5— - x]  -  sin  2j(k-l)xJ  <  e 

8j  (k-ip 


Hence  if 


1  im 


sin[2j(k-l)(l - p— ^ x] 

8j  (k-1 y 


0, 


it  follows  that  x  =  STi/2(k-l),  s  an  integer. 
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We  now  show  that  under  specific  circumstances,  the  tetrachoric  series 

diverges  for  p..  >  p  >  l/(k-l). 

•  J 

Theorem  7-  If  P.-.-  >  P  >  l/(k-l),  T,(x  x  ...,x)  diverges  whenever 

I J  K  »  , 

x  4  sir/2(k-l),  s  an  integer. 

Proof.  Let  Q  (x)  be  the  term  of  T. (x.x, . . . ,x)  with  p. .  =  p,  n . .  =  m. 

»  P  K  1J  1J 

Then 


W*» 


mk(k-l)/2  k  . 
(m!)k(k-l)/2  *  (x)  H(k-l)m-l(x)- 


(3.14) 


We  examine  the  behavior  of  Qm  p(x)  as  m  ®.  From  ErdSl  -i,  Magnus, 
Oberhettinger,  Tricomi  [3],  p.  201, 


"  2n/2 


r(n+1) 


r(§+D 


cos  (  ( 


2n+l 


nf)eX  /4  +  0((^)Zn'h)  ,  (3.15) 


as  n  •+  “.  Replacing  the  gamma  functions  in  (3.15)  by  Stirling's  approximation, 
we  obtain 
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Hnk(X) 


—jr"  cosk((^±l)Jsx  -  -n|)ekx2/4(l+0(n'i2)) 


(3.16) 


Substituting  (3.16)  into  (3.14)  and  again  employing  Stirling's  approximation 


yields 


„  #l  p^(k-l)/22k/2e-kx^/4[(k_1)m.1]k((k-l)m-l)/2 
VpW  '  k( k-1  )/2e-k/2( 2n)  k( k+1  )/4 


.  cosk(  (  k_:^_)_m- 1  ^  x  .  ( ( k-1  jm-l  h)  ( l40(m-^ )  .  (3.17) 


Writing 


((k-l)m-l)k(k‘1)m/2  =  [(k-l)m]k(k'1)m/2(l  -  (-k_V5p)k(k"1),l,/2 


=  ((k-l)m)k(k_1)m/2e"k/2(l+0(m"1)) 


and  substituting  this  into  (3.17),  we  get 


n  (v)  -  ((k-l)p)mk(k'1)/22k/2e~kx  /4[(k-1)fn-l]~k/2 
VpUJ  "  mk<k-l )/4(2")k(k+U/4 


cosk[(--2(k~])m~1)H  x  -  (ii-1Im--1.)7r)(140(m->5)).  (3.18) 

Then,  if  x  f  sir/2(k-l),  s  an  integer,  it  follows  from  Lemma  3  that  there 
exists  a  subsequence  {m^}  such  that 

lim  |Qm  (x)  1  =  oo, 
v-*°  v,p 

whenever  p  >  l/(k-l).  Clearly,  this  precludes  the  convergence  of  the 
tetrachoric  series  for  p. .  =  p  >  l/(k-l),  1  £  i  <  j  <  k  and  consequently 

*  J 

for  Pij  >  p,  1  <  i  <  j  <  k. 

P-ONfi  imtNG  REMARKS.  In  this  paper,  we  have  shown  that  the  tetrachoric 
series  need  not  converge.  We  have  exhibited  some  instances  of  convergence 
and  some  of  divergence.  The  results  of  this  paper  do  not  exhaust  all  the 
possibilities,  that  is,  not  all  the  possible  choices  of  the  p... 

*  J 

1  £  i  <  j  1  k  or  h.,  i  =  l,2,...,k  have  been  treated.  Further  investigation 
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will  be  needed  to  resolve  the  remaining  cases.  However,  we  have  shown 


that  the  tetrachoric  series  should  be  employed  in  applications  only  under 
highly  restrictive  conditions. 
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